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Abstract 

We present an all-electron study of the dynamical density-response function 
of hexagonal close-packed transition metals Sc and Ti. We elucidate various 
aspects of the interplay between the crystal structure and the electron dynam- 
ics by investigating the loss function, and the associated dielectric function, for 
wave-vector transfers perpendicular and parallel to the hexagonal plane. As ex- 
pected, but contrary to recent work, we find that the free-electron-like aspects of 
the dynamical response are rather isotropic for small wave vectors. The crystal 
local-field effects are found to have an impact on the plasmon energy for small 
wave vectors, which gives rise to an interplay with the exchange-correlation ef- 
fects built into the many-body kernel. The loss function lineshape shows a 
significant dependence on propagation direction; in particular, for propagation 
on the hexagonal plane the plasmon hybridizes substantially with fine structure 
due to d-electron transitions, and its dispersion curve becomes difficult to es- 
tablish, beyond the small wave vector limit. The response is calculated in the 
framework of time-dependent density functional theory (TDDFT), based on a 
full-potential linearized augmented-plane-wave (LAPW) ground-state, in which 
the exchange-correlation effects are treated in the local-density approximation. 

Electronic excitations in solids can be studied from the knowledge of the dynamical 
density-response function, which is directly related to observables such as inelastic 
scattering cross sections for fast electrons and hard x-rays. 

For many years, theoretical investigations of valence-electron excitations were car- 
ried out on the basis of the free-electron gas (FEG) model (or jellium model), in which 
valence electrons are described by an assembly of interacting electrons embedded in a 
uniform background of positive charge. The physics of the FEG model is completely 
controlled by only one parameter, the valence-electron density, uq or, equivalently, by 
the electron-density parameter r s which is related to n through l/n = 47r/3(r s ao) 3 , 
ao being the Bohr radius. As we shall illustrate below, in real materials the physics of 
the excitations is more complex; indeed, the loss spectra shows qualitative departures 
from the FEG description. 

In the past few years, ab initio calculations of the density-response function of many 
simple [l]-[7] and noble [8]- [10] metals have been reported. Interestingly, in a number 



of cases it has been concluded that remarkable anomalies revealed by experimental 
measurements arise mainly from band-structure effects [11]. Examples include the fine 
structure in the inelastic X-ray scattering data for aluminum [1], the negative plasmon 
energy dispersion of cesium [3], [5], and the positive dispersion of the linewidth of the 
plasmon in potassium [6] . 

In this paper we utilize an all-electron method to compute the density-response 
function of Sc and Ti starting from the knowledge of an LAPW ground state [6], [7], [12]. 
Our presentation centers on a discussion of the dominant feature of the excitation 
spectrum for these hexagonal closed-packed transition metals for relatively small wave 
vectors — the collective, plasmon-like mode, and its interplay with the effects of the 
crystal lattice; in particular, we assess the impact of the so-called crystal local fields, 
and the dependence on propagation direction. We conclude that the crystal local-field 
effects are not negligible, in the small wave vector limit; this brings about an interplay 
with the effects of exchange and correlation, which thus impact the plasmon dispersion 
for small wave vectors, unlike the case of the FEG. Similarly, the dependence on propa- 
gation direction of the response function is rather significant: for wave vector transfers 
along the hexagonal plane, the plasmon hybidizes with the fine structure caused by 
ci-electron transitions and its dispersion curve becomes difficult to establish beyond the 
small wave vector limit. Comparison is made with the recent pseudopotential-based 
calculations of the density-response function of Sc by Schone and Ekardt [13]; the sur- 
prising result reported by these authors that the plasmon-like excitation is not realized 
for small wave vectors for propagation in the hexagonal plane is not supported by our 
calculations. Our results agree with the intuitive expectation that the free-electron-like 
response is not affected drastically by the propagation direction. 



1 Theory 

The density-response function x(r, r'; uj) [14] of an interacting electron system gives, 
within linear-response theory, the electron density induced in the system when the 
latter is exposed to an external potential V ext (r, uj) through the equation 

p md (r,uj) = J dr' X (r,r';u;)y e *'(r>)- (1) 

In the framework of time- dependent density-functional theory (TDDFT) [15], the 
exact density-reponse function x(r, r'; uj) obeys the integral equation [16] 

x(t,t';u) = x S (r,r';uj) + J dr x J dr 2X S (r, n; u) 

x [u(n - r 2 ) + f X c(r u r 2 ; u)] x(r 2 , r'; u), (2) 

where x 5 ( r ) r '; ^) is the single-particle density-response function of the unperturbed 
Kohn-Sham (KS) system [17], v(r— r') is the bare Coulomb interaction, and f xc (r, r'; cu) 
accounts for all dynamical exchange-correlation (xc) effects. 

For a periodic crystal it is convenient to work in Fourier space, so we write 

i BZ 

X(r,r» = E e^^e-^'^XG^q,^), (3) 

11 k G,G' 



with a similar Fourier expansion for % (r, r';o;). The Fourier coefficients of the KS 
response function, XG,G'( k > u )i can t> e wr itt en as 

1 /k',n — /k'+k,n' 



Xg.c (k, n E E E ^ n _ Ek , +Kn , + h ( u + ir]) 
x(0 k ^ n |e- i ( k+G )• r |^ +kX )(0 k , +kX | e i ( k+G ')• r |0 k , n ), (4) 

where Q represents the normalization volume, G and G' are vectors of the reciprocal 
lattice, and the second sum runs over the band structure for each wave vector k' of 
the first Brilloun Zone (BZ). k , n (r) and £k,n are Bloch eigenfunctions and eigenvalues 
of the KS Hamiltonian of ground-state density-functional theory (DFT) [17], and / k , n 
are the occupation numbers, with due account for spin degeneracy. 

In the above Fourier representation, eq. (2) turns into a matrix equation 



XG,G'( k > = X G ,G'( k > <") + E E XG,G»( k > U ) 

G" G'" 

v G , / (k)5 G „ )G >» + /G» iG '"( k ^)] XG'",G'(k,w), (5) 



which we solve numerically. The size of this matrix equation is a measure of the 
importance of the crystal local-field effects, arising from the inhomogeneity of the 
electronic environment in the periodic crystal potential. For later reference, we note 
that ignoring the crystal local-field effects corresponds to solving eq. (5) while keeping 
only the zero-zero element of the KS response matrix. 

Within the first Born approximation the inelastic scattering cross-section for hard 
x-rays and fast electrons corresponding to wave-vector transfer q = k + G, is propor- 
tional to the dynamic-structure factor 

S(q, u) = -2 h n — ^-rlm \e G ] G (k, u)] , (6) 

where Im e^^ik^) is the so-called energy-loss function and is related to the re- 
sponse function through 

e G | G ,(k,cu) = 5 G:G > +v G (k)xG,G'(k,uj). (7) 

The plasmon-energy dispersion is given by the frequencies at which the real part of the 
dielectric function, e(k, a;), is close to zero and the imaginary part is small. Ime _1 (k, a;) 
is therefore a maximum. Throughout this work the plasmon-energy is taken to be the 
energy location of the main peak of the energy-loss function. 



2 Numerical Implementation 

The key ingredient in the calculation of the energy-loss spectrum is the KS response 
function \ S of eq. (4). For this purpose, we need the KS states which we obtain within 
the local-density approximation (LDA) [18], in terms of a variational expansion in the 
LAPW basis [19]: 

^n(r) = ^E^(k+G)^^(r). (8) 

G 



The LAPW functions ?Pt+G W (. r ) are obtained dividing the unit cell into two regions: 
non-overlapping atomic spheres centered at nuclear sites and the interstitial region 
between the spheres. This allows a faithful description of the localized strong oscil- 
lations near the atomic site, where atomic-like functions are used, and the smooth 
behaviour of the interstitial region, where plane waves are employed. Local orbitals 
are also introduced in order to have an accurate description of the semi-core states. 

The use of the symmetry properties of the crystal allows us to work with Bloch 
states that involve momentum transfers in the irreducible Brillouin Zone (IBZ) instead 
of those in the whole zone. If k is a vector of the IBZ, 

0k,n(r) = R -!k, n (r) = e^-W^jR ■ r + r), (9) 

where R(k) runs over the sub-set of the point group of the crystal that generates the 
star of k and r(R) is the partial lattice displacement corresponding to the space-group 
operation {R, r(R)}. 

The evaluation of the matrix elements in eq. (4) is performed with a code that 
runs in parallel and scales linearly with the number of nodes used. With the matrix 
elements at hand we perform the sum over k-points and over the pairs of bands allowed 
by the occupation numbers. We solve equation (5) using two different approximations 
for f xc . In the random-phase approximation (RPA), f xc = 0. In the so-called adiabatic 
local-density approximation (ALDA), 

(10) 

n=n (r) 

where E xc (n) is the exchange-correlation energy of a uniform electron gas of density 
n, and no(r) is the ground-state density. 

We apply the above formalism to hexagonal closed-packed transition metals Sc and 
Ti, whose electronic configurations are [Ar]3<i 1 4s 2 and [Ar]3<i 2 4s 2 , respectively. The 
lattice parameters used are a=2.95 A and c=4.68 A for Ti and a=3.31 A and c=5.28 A 
for Sc [20]. The ground state is obtained with a cut-off parameter R MT K max =8. Inside 
the atomic spheres, the LAPW wave functions are expanded in spherical harmonics 
Y[ m up to Z=10. The response function was evaluated for q wave vectors perpendicular 
and parallel to the hexagonal plane, using 8 x 8 x 16 and 16 x 16 x 8 BZ samplings 
(the corresponding number of points in the IBZ being 90 and 150, respectively) and 
keeping KS states up to 7.5Ry, thereby ensuring convergence in the energy range under 
study. The 3s and 3p states were treated as semi-core states. 

3 Results 

The upper left (right) panel in Fig. 1 shows the dynamical-structure factor obtained 
within ALDA for Ti (Sc) for energies up to the onset of semicore excitations (~ 32.5 
eV in Ti and ~ 27.7 eV in Sc) along the (001) direction. We have inverted eq. (5) 
using up to 9 shells (55 G- vectors) [21]. For clarity, only results for the scalar inversion, 
corresponding to a lxl matrix calculation with the G = G' = element (dashed line), 
and for matrix inversion using 11 (dotted-line), 23 (solid line) and 35 (open circles) 
G- vectors are displayed. It's apparent that convergence is reached for the energy 
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Figure 1: Top panel: impact of local-filed effects on the dynamic structure function of Ti 
and Sc along the hexagonal axis within ALDA. The values of the momentum transfer are 
q= 4/16(001) and q= 5/16(001) respectively, both in units of 2ir/c. The dotted-dashed line 
corresponds to a FEG calculation with r s =1.92 a.u. (Ti) and r s =2.38 a.u. (Sc). Lower 
panel: dielectric function. 

range under study (the physics of the dynamical response above the M-edge will be 
the subject of another publication); the results shown in Figs. 2 and 3 correspond 
to the converged 35 G-vector calculations [22]. S^q, u) displays a broad plasmon-like 
peak, whose position is shifted upwards (between 0.5 eV and 1.0 eV) by the crystal 
local-field effects. These peaks correspond (see lower panels of Fig. 1) to zero values 
of Re e(q, uj) where Im e(q, uj) is small. 

For comparison, the upper panels of Fig. 1 also show (dotted-dashed line) the 
dynamical structure factor calculated within the FEG model for r s =1.92 a.u. (Ti) 
and r s =2.38 a.u. (Sc) for the same wave-vector as the ab initio calculation, r s being 
the electron density equivalent to that of valence electrons, ?>d l As 2 and 3d 2 As 2 , in Ti 
and Sc respectively [23]. Clearly, as anticipated at the outset, the overall loss functions 
of these transition metals bear little resemblance to the FEG model, other than the 
presence of a broad collective mode; again, the essentially overdamped mode which 
defines the leading feature in Fig. 1 corresponds to a dielectric function whose nature 
is quite removed from the simple, clean, Drude-like zero which defines the textbook 
plasmon. Additional peaks around 9 eV in Ti and around 6 eV in Sc arise from 
transitions involving the (i-electrons and will be discussed elsewhere. 

In order to investigate the dependence on the direction of the momentum transfer, 
we show in Fig. 2 the ALDA dynamic structure factor S(q, uj) for Ti and Sc along 
the (001) and (110) directions for q values of the same order of magnitude. Dashed 
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Figure 2: Comparison of the dynamical structure factor for similar values of the momentum 
transfer perpendicular (001) and parallel to the hexagonal plane (110) for Ti (left) and Sc 
(right) within ALDA. 



and dotted lines represent S(q, u) for q= 3/16(001) and q= 4/16(001), respectively, in 
units of 2n/c. The calculation along the (110) direction (solid line) is for q= 1/16(110), 
in units of 27r/a. The curves along the (001) direction are calculated with a 8 x 8 x 16 
sampling over the BZ and the ones in the (110) direction with a 16 x 16 x 8 sampling 
[24]. As the left panel on Fig. 2 illustrates, for both propagation directions the plasmon 
peak in Ti hybridizes with fine structure due to <i-electron transitions (prominent peak 
and shoulder at its left, respectively). This hybridization is more pronounced along 
the hexagonal plane. For this reason, the identification of a plasmon peak for larger 
momentum transfers along this direction becomes very ambiguous. 

In recent work [13], Schone and Ekardt have reported pseudopotential-based cal- 
culations of the dynamic structure factor for Sc. The intriguing aspect of their results 
is that they suggest that the collective mode is not realized for small wave vectors for 
propagation along the hexagonal plane. By contrast, as shown in Fig. 2 the nature of 
the collective-like response in Sc (and also in Ti) is, for small wave vectors, quite simi- 
lar for propagation directions parallel and perpendicular to the hexagonal plane. Our 
result agrees with intuitive expectations; on that basis we submit that it is correct. 

Fig. 3 shows the plasmon-energy dispersion for both directions within RPA (solid 
symbols) and ALDA (empty symbols). Circles and squares correspond to wave- vector 
transfers that are perpendicular and parallel to the hexagonal plane. As expected, the 
ALDA plasmon energy shifts to slightly smaller values compared to that of the RPA, 
due to a less effective electron-electron interaction. 

An interesting feature in both panels of Fig. 3 is that, in the small wave vector 



Ti 



Sc 



20 
19.5 



> 










19 








<D 




<D 


18.5 


i 

fl 




O 




a 


18 






'Eh 





17.5 
17 



T 



T 



T 



•o° 



o 



•o 
• o 
• o 
•o - 



o 



o 



• 


001 RPA 




13 


o 


001 ALDA 






■ 


110 RPA 






□ 


110 ALDA 







J_ 



_L 



J_ 



J_ 




0.2 0.4 0.6 0.8 

o -1 

q(A ) 



12 



0.2 0.4 0.6 0.8 

o -1 

q (A ) 



Figure 3: Plasmon-energy dispersion for Ti (left) and Sc (right) along the 001 direction 
(circles) and along the 110 direction (squares) within RPA (solid symbols) and within ALDA 
(empty symbols). 



limit, RPA and ALDA do not yield the same result, as it occurs in the FEG. Indeed, 
in the case of the FEG, the Coulomb term dominates over the exchange-correlation 
contribution when q — > 0. In the case of the ab initio calculation, the difference 
between RPA and ALDA arises from the inclusion of crystal local-field effects and its 
interplay with the exchange correlation effects: for small wave vectors, q in eq. (6) 
lies in the BZ and hence, the element we consider after solving the matrix equation 
(5) for xg,c is the G = G'= element; thus, in this case k = q. 

Through the inversion process there is a "mixing" between large wave vector argu- 
ments k+G in x S an d the small wave vector q of the plasmon, causing a shift of the 
collective peak to larger energies. On the other hand, as mentioned before, exchange- 
correlation effects cause a downwards shift, smaller than the former one, provoking 
thus the difference [25]. The partial cancellation between both effects leads to a net 
difference between the converged, ALDA result, and the scalar RPA result, for the 
small-q plasmon dispersion curve. 

Ti and Sc show positive plasmon dispersions and disperse faster along the hexagonal 
axis. For |q| values larger than ~ O.sA^ 1 , the assignment of a collective peak in 
the (110) direction becomes a difficult task. By contrast, this assignment is still 
possible along the (001) direction. This fact ratifies, once more, the importance of the 
underlying crystal structure. 

We note that for Ti, the plasmon energy (~ 19 eV) is in good agreement with the 
one measured in reference [26] for q=0.7 A -1 along (001). 



4 Conclusions 



We have presented results of calculations of the dynamic response function of hep 
Sc and Ti based on a full-potential linearized augmented plane wave method. We 
have found that crystal local-field effects have a non-negligible impact on the plasmon 
energy for small wave vectors causing an interplay between the crystal potential and 
the exchange-correlation effects. We have investigated the dependence on propagation 
direction of the dynamical-structure factor and we have shown that, contrary to a 
recent report [13], the plasmon loss is quite isotropic; however, the lineshape of the 
energy loss shows considerable dependence on the direction of the wave vector transfer: 
along the hexagonal plane, the identification of a dispersion relation for the plasmon 
energy beyond the small wave vector limit, becomes more difficult due to hybridization 
with the fine structure originated by ci-electron transitions. 
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